The origin of flux-flow resistance oscillations in Bi 2 Sr 2 CaCu208+j / : 
Fiske steps in a single junction? 
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We propose an alternative explanation to the oscillations of the flux-flow resistance found in 
several previously published experiments with Bi2Sr2CaCu208+j, stacks. It has been argued by the 
previous authors that the period of the oscillations corresponding to the field needed to add one 
vortex per two intrinsic Josephson junctions is associated with a moving triangular lattice of vortices 
(out-of- phase mode), while the period corresponding to one vortex per one junction is due to the 
square lattice (in-phase mode). In contrast, we show that both type of oscillations may occur in a 
single-layer Josephson junction and thus the above interpretation is inconsistent. 

PACS numbers: 74.50.+r,74.40.+k,74.78.Na 



Recently, a lot of attention has been focused on the 
possible use of single crystals of Bi 2 Sr 2 CaCu2 08+y as 
generators of electromagnetic radiation in the THz range. 
Promising experiments in that direction were reported by 
Hirata et al. and shortly thereafter by Kakeya et al. 
and Hatano et al. 0. These experiments all showed 
oscillations in the flux-flow voltage and flux-flow resis- 
tance when a large magnetic field (of order several Tesla) 
was applied parallel to the ab plane in the presence of 
a small bias current in the c-direction (of order a few 
percent of the critical current). The observed flux- flow 
voltage oscillations typically showed two different oscilla- 
tion periods. At the lowest magnetic fields the period was 
AHt = $o/(2sL) corresponding to one extra flux quan- 
tum per two layers in the stack. Here $o is the magnetic 
flux quantum, s and L are the thickness and the length 
of the junction, respectively. At higher magnetic fields 
there was a transition to a period AHs = $o/(sL) , i.e. 
corresponding to an extra flux quantum in every layer. 

These exiting experiments were interpreted both an- 
alytically and numerically 0, H, by several au- 
thors. One of the important questions to be answered 
was whether the flux lattice correspond to a triangular 
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lattice (anti-phase ordering with possible cancellation of 
the sum voltage at the junction end) or a square lat- 
tice (in-phase ordering leading to a large sum voltage 
at the junction end). Figure ^ shows schematically a 
F^S^CaCu^Os+j, stack with flux ordering in a triangu- 
lar lattice and square lattice. Obviously, the latter case 
is highly preferable for applications to THz generation of 
electromagnetic waves; however simple intuition would 
suggest triangular ordering since fluxons of same polarity 
naturally repel each other. An intuitive interpretation of 
the experimentally observed oscillations in the flux-flow 
voltage would suggest that an oscillation period AHt 
corresponds to triangular ordering while a period AHg 
corresponds to a square lattice. Numerical simulations 
0- have shown that both triangular and square lat- 
tices are possible, but their relation to regions of AHt 
oscillations and regions of AHs oscillations is not simple 
and details are still a matter of debate. 



In this paper, we present fairly standard numerical sim- 
ulations corresponding to the simplest case of a single- 
layer Josephson junction. Even for this case, where there 
is no triangular nor square lattice ordering (as our stack 
consists of only one junction), we find that both the AHt 
and AHs periods appear much the same way as in exper- 
iments and numerical simulations for I^S^CaCu^Og+j, 
stacks with many layers and flux lattice ordering. Af- 
ter presenting the numerical simulations we provide a 
qualitative explanation in terms of the well-known Fiske 
modes lg. 



FIG. I: Schematic drawing of the triangular lattice (a) cor- 
responding to the out-of-phase fluxon mode and the square 
lattice (b) associated with the in-phase fluxon mode. 



The system under investigation is described by coupled 
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sine- Gordon equations of the form Q: 
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Here a = (l/i?)[?i,/(2e/oC)] 1/2 is the dissipation param- 
eter (R, Io, and C are the normal resistance, the crit- 
ical current and the capacitance per unit length, re- 
spectively), 7 is the current normalized to the critical 
current Iq of the individual junctions. The normalized 
coupling term among the junctions in the stack reads 
S = -A L /Ksinh(t/A L )], where d! = d + 2A L coth(t/A L ), 
t is the thickness of one superconducting layer, and Al 
is the London penetration depth Time t is nor- 

malized to the inverse of the Josephson plasma fre- 
quency uj p = [2eIo/(HCj )] 1 ^ 2 and spatial coordinate 
x is normalized with respect to the Josephson length 
Aj = [?i/(2e/io^o)] 1 ^ 2 - The magnetic field gives rise to 
boundary conditions of the form 
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Index i = 1, . . . , N stands here for the junction number 
in the stack and £ is the normalized length of the system. 

Figurc[21shows the result of a simulation of the current- 
voltage characteristics for only one junction in the stack, 
i.e. N = 1. It shows the flux- flow branch of the junction 
with normalized length £ — 40 and damping constant a = 
0.1 placed in magnetic field h = 4. The characteristics is 
calculated by rising the bias current 7 from zero to 0.85 
and then decreasing it back to zero. The voltage V is 
given in normalized units chosen such that the voltage 
spacing between neighboring Fiske steps 



AV = 



2L 



(4) 



is equal to unity (here c is the Swihart velocity and 
L = £Xj is the physical length of the junction). The 
current- voltage characteristics displays fine structure due 
to the Fiske steps, which gather around the flux-flow volt- 
age, also known as Eck peak ^(|- The hysteresis in this 
voltage region is due to the coexistence of several Fiske 
resonances at a given current 7. Some parts of these Fiske 
resonances are located inside the hysteresis and are not 
displayed on the plot. 

Using the above junction parameters, we calculated 
the dependence of the flux-flow resistance on magnetic 
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FIG. 2: Current-voltage characteristics of a long junction 
(N = 1) with parameters £ — 40, a = 0.1, and h — 4. Arrows 
indicate switching between branches for rising and decreasing 
bias current 7. 




magnetic field h 

FIG. 3: Differential flux-flow resistance dV/d7 versus mag- 
netic field h for the single junction (N = 1) with I = 40 and 
a — 0.1 at constant bias current 7 = 0.3. 



field. Figure presents the differential flux-flow resis- 
tance dV/d'y versus magnetic field h at the fixed bias 
current 7 = 0.3. It clearly shows flux-flow resistance os- 
cillations. In order to relate the period of oscillations 
with the number of vortices in the junction we show a 
grid in magnetic field with a period Ah = 2tt/L£, which 
approximately corresponds to adding one vortex in the 
junction. This can be seen from the simple fact that 
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magnetic field h 

FIG. 4: Flux- flow resistance V/7 (thick curve) and AV/A'y 
(thin curve) versus magnetic field h for the same junction as 
in Fig. |3]at a lower bias current 7 = 0.2. 

at the critical field h = 2 the normalized spacing be- 
tween vortices penetrated into the junction is equal to 
7r, and their number rises proportionally to h. In Fig. El 
we see that most oscillations have a characteristic period 
in h corresponding to adding a half flux quantum into 
the junction, with a tendency of doubling the period at 
higher fields. 

As another illustration, in Fig. 0] we show both static 
resistance V/7 and differential resistance dV/d7 versus 
magnetic field h for the same junction but at a lower 
bias current 7 = 0.2. The oscillations of dV/d7 (thin 
curve) have two characteristic periods, which are found 
at different ranges of magnetic field h. The oscillation pe- 
riod at low fields corresponds to about or less than half 
flux quantum, while at higher fields we find very clear 
oscillations which account for one flux quantum into the 
junction. The crossover from one regime to another oc- 
curs at magnetic field h w 3. Below this field the differ- 
ential resistance at 7 = 0.2 is lower as it is determined 
by the Eck peak (see also Fig. |2J) composed of individ- 
ual Fiske steps. At h > 3 the Eck peak shifts to higher- 



currents and the differential resistance levels at the re- 
sistive slope determined by the loss parameter a of the 
junction. The static junction resistance (upper curve in 
Fig. also changes in this range but its oscillations are 
much less pronounced and can only be clearly seen on 
the magnified scale. 

We suggest the following explanation for the two os- 
cillation periods of the flux-flow resistance in a single- 
barrier Josephson junction. At high enough magnetic 
field the resistance at a low bias current follows oscil- 
lations of the critical current of the junction and thus 
have a characteristic period of one flux quantum. At the 
same time, the resistance measured at a higher current 
(or lower field but the same bias current) follows the os- 
cillations due to the Fiske steps, which envelope is associ- 
ated with the Eck peak (sometimes called as flux-flow or 
velocity matching step). The matter is that - for a single 
junction - the Fiske steps induce a variation of the re- 
sistance with characteristic period corresponding to half 
flux quantum. 

Fiske modes in a Josephson junction 

8] 

are linear cav- 
ity type excitations with resonance angular frequencies 
given by 

2n 

u n =n—, n = 1,2,3,... (5) 

in normalized units. The corresponding to wave vectors 
are 

fen = n— . (6) 

In experiments these Fiske modes are visible as current 
singularities in the current-voltage curve with a voltage 
spacing given by Eq. (@J. The amplitude of the Fiske 
steps oscillate with the magnetic field in a typical Bessel 
function like pattern such that the even numbered steps 
have maxima together with maxima in the critical cur- 
rent, while maxima in the odd numbered steps corre- 
spond to minima in the critical current. Refs. 0, 0] 
give an approximate analytical form for the current- 
voltage curve for a single Josephson junction with Fiske 
steps. The current-voltage characteristics is approxi- 
mately written as 0] : 
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This approximation neglects high-order nonlinearities in well the case of short junction, i.e. t < 1. However, the 
the junction and is originally expected to describe comparison with long junction data made by Cirillo et 
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FIG. 5: Differential flux- flow resistance dV/d7 at constant 
bias current 7 = 0.2 versus magnetic field h calculated directly 
from Eq. Q with parameters I = 40 and a = 0.1. The voltage 
is multiplied by a factor £/n to make its scale identical to the 
voltage normalization used in the previous figures. 

al. [13 shows that Eq. Q can approximately account 
for the shape and for the maximum current modulation 
of the Fiske singularities in long (£ 3> 1) junctions when 
the field penetration overcomes Meissner shielding, i.e. 
at h > 2. 

The first term in Eq. Q) represents the Ohmic part of 
the current- voltage characteristics, while the second term 
gives an infinite series of equidistant resonances. The 
height of the resonances is modulated by a slowly varying 
amplitude factor and a fast Fraunhofer amplitude factor 
|l2| . The Fraunhofer factor emphasizes the resonance 
closest to lu — h and drops off fast: If hi is an even 
multiple of 7r, the odd numbered Fiske steps are enhanced 
and if hi is an odd multiple of 7r, the even numbered Fiske 
steps survive. 

Equation JJJ gives the current-voltage curve contain- 
ing Fiske steps with the magnetic field as a parameter. 
If we instead assume a fixed bias current 7 and vary the 
magnetic field, Eq. JJJ) expresses the voltage V oc ui os- 
cillations in an implicit form. In order to compare this 
analytical form with our simulations, we solved Eq. (JZJ) 
for w at a given 7 numerically. The obtained the depen- 
dence of the flux-flow resistance dV/d'y at constant bias 
current 7 = 0.2 versus magnetic field h is presented in 
Fig. The qualitative agreement between Figs. EH and 
0] obtained by full numerical simulation of the perturbed 
sine-Gordon equation and Fig. [3] emerging from the an- 
alytical formula Q) is strikingly good. Figure [5] clearly 



displays two characteristic periods of oscillations, namely 
half flux quantum oscillations at low fields and one flux 
quantum oscillations which become very explicit at high 
fields. The intermediate field range shows a complicated 
beating between two periods. 

Thus, odd and even numbered Fiske resonances in 
the last term in Eq. (JTJl produce the "magic" half-flux- 
quantum oscillations corresponding to the magnetic field 
period AHt = $0/ (2sL) even in a single Josephson junc- 
tion. Although we investigated here Fiske steps with 
N = 1 , we note that Fiske steps are also present in stacks 
with N > 1 [l^. Thus we are lead to suggest that also 
for N > 1 the flux-flow voltage oscillations have their 
origin in the Fiske mode excitations. 

We conclude that for Bi2Sr2CaCu20s+ a stacks the 
flux-flow voltage oscillations with two different periods 
in a magnetic field have their origin in the Fiske mode 
excitations. Thus the flux lattice ordering in either tri- 
angular or square lattice is not directly related to the 
two periods of the oscillations. We speculate that Fiske 
modes also existing in stacks indirectly play a role for the 
flux lattice formation. 

We would like to acknowledge discussions with T. 
Hatano, S. Kim, A. E. Koshelev, S. Madsen, M. R. 
Samuelsen, H. B. Wang, and T. Yamashita. 



[1] S. Ooi, T. Mochiku, and K. Hirata, Phys. Rev. Lett. 89, 
247002 (2002) 

[2] I. Kakeya, M. Iwase, T. Yamamoto, and K. Kadowaki, 
cond-mat/0503498 (2005); I. Kakeya et al. Physica C 
378,406 (2002) 

[3] T. Hatano et al., IEEE Trans Appl. Superconduc, to be 
published 

[4] A.E. Koshelev, Phys. Rev B 66, 224514 (2002) 

[5] M. Machida, Phys. Rev. Lett. 90, 037001 (2003) 

[6] M. Tachiki , M. Iizuka, K. Mi nami, S. Tejima, and H. 

Nakamura, |cond-m at/0407052 (2004) 
[7] N.F. Pedersen and S. Madsen, IEEE Trans Appl. Super- 
conduc, to be published 
[8] M. D. Fiske, Rev. Mod. Phys. 36, 221 (1964) 
[9] S. Sakai, P. Bodin, and N.F. Pedersen, J. Appl. Phys. 

73, 2411 (1993) 
[10] R. E. Eck, D. J. Scalapino and B. N. Taylor, Phys. Rev. 

Lett. 13, 15 (1964). 
[11] I. O. Kulik, Pis'ma Zh. Eksp. Teor. Fiz. 2, 134 (1965) 

[Sov. Phys. JETP Lett. 2, 84 (1965)]; I. O. Kulik, Zh. 

Tekhn. Fiz. 37 157 (1967) [Sov. Tech. Phys. 12, 111 

(1967)] 

[12] M. Cirillo, N. Gr0nbech- Jensen, M. R. Samuelsen, M. 

Salerno, and G. V. Rinati, Phys. Rev. B 58, 12377 (1998) 
[13] S. Kim, S. Urayama, H. B. Wang, T. Hatano, M. Na- 

gao, S. Kawakami, Y. Takano, and T. Yamashita. Poster 

Presentation in the conference Rlasma-2004' Tsukuba, 

Japan, November 2004 



